Superconductivity in Inhomogeneous Hubbard Models 
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We present a controlled perturbative approach to the low temperature phase diagram of highly 
inhomogeneous Hubbard models in the limit of small coupling, t' , between clusters. We apply 
this to the dimerized and checkerboard models. The dimerized model is found to behave like a 
doped semiconductor, with a Fermi-liquid groundstate with parameters (e.g. the effective mass) 
which are smooth functions of the Hubbard interaction, U. By contrast, the checkerboard model 
has a nodeless d-wave superconducting state (preformed pair condensate, d-BEC) for < U < U c , 
which smoothly crosses over to an intermediate BCS-like superconducting phase (d-BCS), also with 
no nodal quasi-particles, for \U — U c \ < 0(t'), which gives way to a Fermi liquid phase at large 
U >U C = 4.58. 

PACS numbers: 74.20.Mn, 71.10.Fd 



In this paper, we report a study of inhomogeneous 
Hubbard models, Eq. (JTJ) , in which the lattice is broken 
up into a periodic array of weakly coupled clusters. We 
focus on the case of small "doping," |x| <§C 1, where l±x is 
the number of electrons per site, and on the limit in which 
the coupling between clusters, t', is much less than the 
relevant energy scales within a cluster. Exploiting these 
small parameters, we obtain a well controlled solution of 
the ground-state and low energy excited states. 

There are two purposes of this study: 1) With t' as 
the small parameter, we can trace the non-perturbative 
evolution of the electronic structure as a function of the 
strength of the Hubbard interaction, U, all the way from 
the weak to the strong coupling limit. 2) In light of the in- 

■ creasing evidence^ that some form of self-organized elec- 
| tronic inhomogeneity is widespread in the cuprate super- 
. conductors, it is reasonable to explore the circumstances 

■ in which high temperature superconductivity from purely 
repulsive interactions may be enhanced by certain forms 
of inhomogeneity2i£i^£i£i 

We develop a general strategy for such problems which 
we apply explicitly to the case of the dimerized Hubbard 
model (Fig. 1(a)) and the checkerboard Hubbard model 
(Fig. 1(b)). In both these cases, the undoped (x = 0) 
"parent" Mott insulating system has a unique, insulating 
ground-state with a large spin-gap, A s : 

1) The doped dimerized Hubbard model has a spectac- 
ularly featureless phase diagram. At energies small com- 

' pared to A s , it behaves like a doped semiconductor, with 
, a small Fermi surface enclosing a Luttinger volume equal 
to x, and with an effective mass which changes by a fac- 
tor of 2 as U is increased from U = to U 1. If some 
form of attractive interaction is added to the dimerized 
Hubbard model, such as an additional nearest-neighbor 
exchange energy, J, there is a transition to a singlet su- 
perconducting phase for sufficiently large J, as indicated 
in Fig. 2. However, the superconducting state has mixed 
d- and s-wave symmetry, and a full gap to quasiparticle 
excitations. 

2) The doped checkerboard Hubbard model exhibits four 
distinct zero-temperature phases as a function of U, as 
shown in Fig. 3(a): For < U < U c » 4.58, the 
system is superconducting while for U > U c it is a 
non-superconducting Fermi liquid. The superconducting 



FIG. 1: Schematic representation of the (a) dimerized and (b) 
checkerboard models. The hopping matrix elements are t = 1 on 
the bold lines and t' <C 1 on the thin lines. 

state has d-wave symmetry. Despite this, the quasiparti- 
cle spectrum in the superconducting state is fully gapped, 
even in the range of U within 0(t') of U c , where BCS- 
like d-wave superconducting state occurs (explain below) . 
The Fermi liquid phase is unusual in the sense that there 
are two degenerate bands ("flavors") of fermions with 
plus and minus chirality (in addition to the two spin po- 
larizations). At very large U > 18.6 there is an additional 
transition to a Fermi liquid with spin 3/2 quasiparticles. 

It is also worth noting that inhomogeneous systems 
naturally exhibit precursor superconducting correlations 
well above the actual superconducting T Cl reminiscent of 
some of the pseudo-gap phenomena seen in underdoped 
cuprate superconductors 8 . This tendency is apparent in 
our results, where under many circumstances, the pair- 
ing scale is determined by interactions within a cluster, 
while T c is proportional to (t 1 ) 2 . As a consequence, pair- 
ing persists to a temperature T pa i r ~ (t 1 ) , while the 
superfluid density, and hence the phase ordering temper- 
ature is parametrically smaller, T c ~ (t') 2 , as can be seen 
in Fig. 4. 



I. THE INHOMOGENEOUS HUBBARD MODEL 

While ideally we would like to consider a system in 
which any inhomogeneity is self-organized, in the present 
paper the inhomogeneity is introduced explicitly from the 
beginning. We therefore consider the Hubbard model on 
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a square lattice 
H = - £ U 
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where (r, r') indicates nearest-neighbor sites, and 
creates an electron on site r with spin polarization a = 
±1. The term with U > represents the on-site repul- 
sion between electrons and n Tta — c\. a c r ^. The usual 
(homogeneous) limit of this model is obtained by taking, 
t r y = t = 1, where the final equality defines our units of 
energy. In the inhomogencous versions of this model we 
consider, the lattice is broken up into a set of periodically 
repeated disconnected clusters, with t Ttr i — 1 for nearest- 
neighbor sites within a single cluster, and t T y — t' <C 1 
for nearest-neighbor sites belonging to distinct clusters. 
Even the inhomogeneous version of this model is particle- 
hole symmetric; we will discuss the case of a concentra- 
tion x of doped holes, but the same results apply for the 
same concentration of doped electrons. 

To zeroth order in t' , the Hamiltonian can be solved 
for arbitrary U by diagonalizing it on a single cluster. We 
then use low order (near) degenerate perturbation theory 
(in powers oft') to derive an effective Hamiltonian, H e *\ 
which operates in the reduced Hilbcrt space spanned by 
the direct products of the low-energy eigenstates of the 
isolated cluster. This is a standard procedure, precisely 
analogous to that used to derive the t-J model from the 
large U limit of the Hubbard model^. For all the clusters 
we consider here, the groundstate of the isolated undoped 
cluster (with one electron per site), is a spin singlet with 
a finite spin gap A s . 

For small x, most clusters must still be in their ground- 
state, so H e *' operates in a very much smaller Hilbcrt 
space than the starting space. Moreover, defining the 
unique ground-state with one electron per site to be the 
vacuum state of iJ e ", it is clear that it can typically be 
recast as the Hamiltonian of a dilute gas of excitations. 
This is the key feature that makes the problem tractable 
in the stated limit of small x and small t' . 

To construct the low energy Hilbert space, we need 
to compute the spectrum of an isolated cluster with dif- 
ferent numbers of doped holes. The eigenstates of each 
cluster can be identified by their symmetry related quan- 
tum numbers: the number of doped holes, Q (Q = 
refers to the case of one electron per site) the total spin, 
S, and those related to the spatial symmetries. For the 
dimer, the states are odd or even under reflection. The 
isolated square has the same four-fold rotational symme- 
try, C4, as the uniform lattice so the states can be labeled 
by spectroscopic labels "s" (even under rotation by n/2) 
"d" (odd under rotation by n/2) and "p x ±ip y v (changes 
phase by ±7r/2 under rotation by ir/2). In each charge 
sector, so long as there is a "large" (order 1) gap, the ex- 
cited states can be safely eliminated from the low energy 
sector. Where there is a level crossing within the isolated 
cluster, we need to be a bit more careful. 

Isolated dimer: For the isolated dimer with Q = 
or Q = 2, there is a unique S — 0, even parity ground 
state separated by a large gap from the first excited state. 
The Q = 1 ground-state is a S = 1/2 even parity doublet 
again with a large gap. 
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FIG. 2: (a) Zero-temperature phase diagram of the dimerized Hub- 
bard model for small x, and (b) phase diagram of H e ^ on the 
dimerized lattice as functions of A p /r and g/r. Two phases, the 
Fermi liquid phase (FL) and the singlete superconducting phase 
(s+d SC), can be obtained in the effective thoery. However, as 
shown by the dotted curve with arrows in (b), the trajectory of the 
Hubbard model with fixed t' /t (=0.005) as a function of increasing 
U only passes through FL. As a result, only FL can be seen in (a). 

Isolated square: For the isolated square with Q = 2, 
there is a unique S — ground state with s-wave sym- 
metry separated by a large gap from the first excited 
state. For Q = and with U = 0, there is a large 
(6-fold) ground-state degeneracy. This degeneracy is 
lifted* at non-zero U, and the resulting ground-state is an 
S = singlet with c?-wave symmetry However, for small 
enough U the gap to the excited states is small - the 
splitting between the lowest lying singlet and triplet state 
(the "spin-gap") is 0(U 2 ) for small U . In the present pa- 
per, when dealing with the checkerboard lattice, we will 
assume U \ft/ so that the gap can be treated as "big" , 
but the small U limit is probably worth revisiting in the 
future. As pointed out by Trugman and ScalapinoiS, an 
important consequence of the distinct spatial symmetries 
of the Q = and Q — 2 ground states is that the pair cre- 
ation operator that connects these two-states has e?-wave 
symmetry. The Q = 1 spectrum of the isolated square 
is a bit more complex: For U < Ut — 18.6, the ground- 
state is a spin and orbital doublet, with S — 1/2 and 
p x ± ip v symmetry. This has the consequence that the 
quasiparticles carry an orbital "flavor" index, A = ±1, in 
addition to the usual spin polarization index, a = ±1/2. 
For U > Ut, the ground state, in accordance with Na- 
gaoka's theorem^, is a S = 3/2 s-wave state. Except in 
the vicinity of U — Ut (where the gap is 0(\U — Ut\), 
the gap to excited states is again large. 

Other than the stated level crossings, the precise de- 
pendence of the energies of the various states is not im- 
portant for present purposes. The energy of the Q = 
groundstate can be absorbed into an overall constant con- 
tribution to the effective Hamiltonian, Eq, and the energy 
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of the Q = 1 state into a redefinition of the chemical po- 
tential. There is one important combination of energies, 

A p = 2E(1) - E{2) - E(0), (2) 

where E(Q) is the ground-state energy for given Q. This 
has the interpretation of the pair binding energy: a pos- 
itive A p signifies an effective attraction between doped 
electrons or holes in the sense that for two doped holes, 
it is energetically preferable to place both on one cluster 
than to place one on each of two clusters. For the isolated 
square, A p is positive (pair binding) for U < U c « 4.58, 
and negative for U > U c . For the isolated Hubbard 
dimer, A p < for all [/, and indeed it vanishes (lin- 
early) only at U = 0. This is the reason the dimerized 
Hubbard model does not superconduct. 



II. THE EFFECTIVE BOSON-FERMION 
MODEL 

It is now straightforward to obtain the effective Hamil- 
tonian on the cluster lattice to first order in t' taking the 
unique ground-state of the undoped system as the vac- 
uum state (See appendix for the derivation): 

H eff = E a +Y^ {-\n b] - (i[2n bj +n fj }) 
j 

- EE (".•<v\ '',~.,. A ",'.-'.A • ii-<-) 

(ij) <r,XX' 

+ E (pxx'^Hjbt [ai,-[,xaj,i,x> - (?<->!)] + H.c)j 

m,xx' 

+ r^fnj.+ns.X^ + riM-l) (3) 

i 

where a\ ■> creates an one-hole fermion on the zth cluster 
with spin polarization a and (for the checkerboard case) 
flavor index A and b\ creates a hole-pair boson on the 

ith cluster; n.&j = bjbj and n/j = J2 a \ a \ a A a j>>A are i 
respectively, the boson and fermion densities on cluster 
j. The coupling constants, r, g (both proportional to 
t'), and our old friend A p represent the effective hop- 
ping of one- hole fermions, the fcrmion-boson Andreev 
coupling, and the energy difference between one boson 
and two fermions respectively, (ij) represents a pair of 
nearest neighbor clusters, and 4>ij is a geometric factor 
discussed below. The effective Hamiltonian operates in a 
constrained Hilbert space where + ntj = or 1, but 
equivalently, it can operate in an unconstrained Hilbert 
space with the constraint imposed dynamically by tak- 
ing the limit U°° — ► oo. The chemical potential is, of 
course, an implicit function of the doping concentration 
x of the original model, obtained by inverting the rela- 
tion '^2 i (nfi + 2nbi) = Mx, where M denotes the number 
of lattice sites in the original Hubbard model. 

This effective Hamiltonian is not only built on a highly 
reduced Hilbert space, but also has common structure for 
both models of our current study. It has a fermion band 
(or multi-bands with index A), a boson band (with in- 
finite bare effective mass) and the interactions between 
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FIG. 3: (a) Zero-temperature phase diagram of the checkerboard 
Hubbard model for small £(=0.025) and (b) of H e ' ' on the checker- 
board lattice. The two superconducting phases both have d-wave 
symmetry and no nodal quasi-particles, but one is BCS-like (d- 
BCS) and the other is a Bose condensate of preformed two-particle 
bound-states (d-BEC). The dashed lines indicate a crossover rather 
than a sharp phase transition. Notice that the boundary (not 
drawn here) for the presence of a two-particle bound state in Eq . J3t 
is close to the crossover curve. The two Fermi liquid phases (FL 
and FL') are distinguished by the quantum numbers of the quasi- 
particles. The dotted curve with arrows in (b) represents the ef- 
fective trajectory corresponding to the Hubbard model with fixed 
t'/t (=0.005) as a function of increasing U. 

them which converts a pair of fermions to a boson and 
vice versa. The effective theory is analogous in form 
to the so-called "Boson-Fermion" model which has been 
studied by several people 12 . However, because in the 
present case, this model is derived as the low-energy 
effective field theory from the inhomogeneous Hubbard 
model, we are lead to study it in particular limits (espe- 
cially i€ 1) that were not the focus of previous studies. 

For the dimerized model, the dimers explicitly break 
the C4 symmetry of the underlying lattice, and hence 
H e ff only has C2 symmetry. Explicit evalution of the 
first order perturbation theory leads to the expressions 
A p = -2t{l + t&nO), g = (cos8 + s'm9)t' , and for 
nearest neighbor dimers in the same column of dimers, 
Tij = (1 — sin2f?)£'/2 and 4>ij = 1) while for dimers in 
neighboring columns, ry = (1 — sin20)t'/4 and <f>ij = 1/2 
where tan26> = —At/U. 

The checkerboard model preserves the C4 symmetry 
of the underlying lattice. The eigenvalue problem on the 
isolated square is more complex, so the expressions for 
the effective couplings are somewhat complicated. Some 
details are shown in the appendix. To give representative 
values, we give expressions for U in the neighborhood of 
U = U c where Tij t xx' = iT~4>ij£\\i with t = [ao + tt\(U — 
U c ) + . . ]t' for all pairs of nearest-neighbors (an = 0.2217, 
ai = -0.0081), and gw = g[l - 5\,\>] where g = \J3 + 
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(3i(U -U c ) + .. ]t' and (3 = 0.3123, fa = -0.0111, and 
as a consequence of the d-wave symmetry of the pair 
creation operator 4>ij = 1(— 1) for neighboring squares 
along the x-axis (y-axis). (Here e\x is the Levi-Civita 
tensor.) 



III. T = PHASE DIAGRAMS OF THE 
DIMERIZED AND THE CHECKERBOARD 
HUBBARD MODELS 

Let us now consider the phase diagram of H e ff. As 
mentioned above, for small x, the density of excitations 
is small, so that two particle collisions are rare (order 
x 2 ) and multiparticle processes even rarer. Thus, when 
|A p | 3> t', the phases are pretty obvious: 

1) For A p < 0, the bosons appear only as virtual states, 
leading to a weak induced attraction between fermions 
of order V md ~ g 2 /\A p \. However, because of the hard- 
core repulsion between fermions, the net interaction is 
repulsive, and hence there is a Fermi liquid phase. (There 
could be some form of Kohn-Luttinger instability^ at 
exponentially low temperature, but we will not worry 
about this.) For the dimerized model, the Fermi surface 
is a small ellipse closed about k = (7r/2a,7r/a) where a 
is the lattice constant of the original square lattice (FL 
in Fig. 2(a)). In the checkerboard model with U c < U < 
Ut, there are two small Fermi surface circles (due to the 
flavor degeneracy); one is closed about k = (0, 7r/2a) and 
the other about k = (7r/2a,0) (FL in Fig. 3(a)). For 
U > Ut, the system forms a Fermi liquid of spin 3/2 
fermions (FL' in Fig. 3(a)). There is a (presumably first 
order) transition between these two sorts of Fermi liquid 
which occurs at U = U T + 0{t'). 

2) For A p > 0, it is the charge Q = 1 quasiparti- 
cles that are virtual excitations. Integrating them out 
generates an effective nearest-neighbor hopping matrix 
T b^ ~ 5 2 /A p , and a nearest-neighbor boson-boson re- 
pulsion of the same order—. Thus, in this limit, the sys- 
tem has a singlet superconducting groundstate. In the 
dimerized model there is no other symmetry classifica- 
tion of the superconducting state possible (it is an admix- 
ture of d-wave and s-wave states), but for the checker- 
board model, the superconducting state inherits the d- 
wave symmetry of the constituent bosons. This d-wave 
symmetry could be observed in any of the phase sensitive 
measurements that have been discussed in the context of 
the cuprates themselves. However, it is important to real- 
ize that, for both the dimerized and checkerboard model, 
there are no gapless spin 1/2 excitations, and hence no 
nodal quasiparticles. In the language of BCS theory, one 
can think of the d-wave state in the checkerboard lattice 
as being in a strong-coupling limit in which the chemical 
potential has passed below the band-bottom. 

There are several features of the thermal evolution of 
the system in this limit that warrant mention. The first 
is that, since it is a system of preformed bosons, T c 
is determined by the zero temperature superfluid den- 
sity, T c ~ t^^x ~ (i') 2 A~ 1 a:. Of course, since the 
model is two dimensional, this also means that there 
should be a vortex gas regime above T c with a Kosterlitz- 
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FIG. 4: Schematic finite temperature phase diagram of the 
checkerboard model. Notice that this figure is not drawn to scale 
and the x-axis has been offset from zero. As in Fig. 3(a), there 
is a crossover within the superconducting phase from a a!-BEC to 
a d-BCS limit. For large enough U there is a FL phase. Above 
T* ~ 0(t' ), indicated by the top dashed curve, the effective theory 
breaks down. The lower dashed curve, T p , signifies the crossover 
between a normal Bose liquid of preformed pairs for T < T p and a 
high temperature state with little electron pairing for T > T p . The 
range T c < T < T p is a pseudo-gap regime. This can be a broad 
regime since the energy scale of coherence between preformed pairs 
is 0(t' 2 ) as indicated by the solid curve. 

Thouless transition to a phase with quasi-long-range su- 
perconducting order. In addition, there can be two other 
pseudo-gap scales apparent up to temperatures which are 
parametrically larger than T c : Firstly, T p ~ A p /| m(2x)| 
is the characteristic temperature at which pairs ther- 
mally dissociate. Above this, there is a temperature 
T* ~ 0(t' ), at which excitations beyond those in the ef- 
fective model become significant. Our finite temperature 
results discussed above are summarized in a schematic 
diagram, Fig. 4. 

3) The only part of the phase diagram which is at all sub- 
tle is where A p ~ t', where both bosonic and fermionic 
excitations participate in the low energy physics. How- 
ever, this is also in some ways the most interesting 
regime. In the first place, approaching this regime from 
Ap > 0, we see that the superconducting T c oc 1/A P 
reaches its maximum in this regime. Moreover, it is 
only in this regime that we could possibly have a super- 
conducting state with gapless nodal quasi-particles (al- 
though none occur in this model, unfortunately). 

The feature that makes this regime tractable is that 
the transition from the Fermi liquid to the superconduct- 
ing phase is BCS-like. This is because, for small x, it 
is only the pairwise interactions between quasi-particles 
that are relevant. The effective interaction is the sum 
of the induced attraction, V , and the hard-core re- 
pulsion between quasiparticles, which, treated in the T- 
matrix approximation, gives rise to an effective repulsion 
Vo ~ t'. Where the sum, Vt of these two terms switches 
from being a net repulsive to a net attractive interaction, 
it passes through a transition point at which Vt = 0. 
Even though the Fermi energy is only of order t , ar- 
bitrarily near the transition point, \Vt\ Ep, justify- 
ing the use of BCS mean-field theory. For H e * * on the 
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dimerized lattice, since the superconducting state has a 
substantial s-wave component, we generally expect the 
superconducting state near the critical value of A p to be 
nodeless, and thus adiabatically connected to the super- 
conducting phase at large A p . On the other hand, for 
the checkerboard lattice, the superconducting state has 
d-wave symmetry, but a feature of the Fermi surfaces en- 
closing (ir/2a, 0) and (0, 7r/2a) of the BZ is that the line 
of gap zeros does not intersect them, and consequently 
there are no gapless quasi-particles in this phase. As A p 
gets larger, the superconducting state smoothly evolves 
from a BCS-like region, where quasi-particle fluctuations 
are dominant (d-BCS), to a regime of pre-formed pairs, 
where phase fluctuations of the order parameter domi- 
nate the physics (d-BEC). This is quite distinct from a 
rf-BCS phase with nodal quasi-particles. 

To obtain the explicit phase diagram for H & '* pre- 
sented in Fig. 2 and 3, we have made use of a mean- 
field approximation which, none-the-less, reproduces all 
the qualitative features required by the above qualita- 
tive analysis. Specifically, we perform a standard mean- 
field decomposition of H e ^ , and solve self-consistently 
for the anomalous expectation values (bp and (at ^a]^). 
To avoid repetition, we would like to refer the readers 
to reference o 12 i 14 . In this analysis, the hard-core con- 
straint for the bosons is treated exactly, but the onsite 
repulsion involving the fermions is treated as an effective 
repulsion, Vq ~ t'\ predictably, the results do not depend 
substantially on Vq. 

On the dimer lattice, the phase diagram of H e '' as a 
function of A p /r and g/r, shown in Fig. 2(b), exhibits a 
Fermi liquid and a nodeless d+ s superconducting phase. 
If we consider only the allowable region of parameters 
that can be derived from the dimerized Hubbard model, 
these always lie within the Fermi liquid phase as indicated 
by the dotted curve in Fig. 2(b). The arrow on the curve 
represents the direction of change as U increases from to 
a large value. If we had instead considered the dimerized 
t-J model, then for J/t > 2 - 22.9(i'/t) the system is 
superconducting. 

The phase diagram of H e ^ on the checkerboard lattice 
is shown in Fig. 3(b). The dotted curve in Fig. 3(b) is 
the trajectory through the phase diagram derived from 
the checkerboard Hubbard model with fixed t' = 0.005, 
where the arrow indicates the direction of change as U 
increases from to a large value not too close to Ut- It is 
worth emphasizing that, in this case, superconductivity 
with d-wave symmetry is obtained from purely repulsive 
interaction and can survive in the region with slightly 
negative A p (~ 0(-f)). 



IV. PERSPECTIVE 

In the present study, the translational symmetry of 
the lattice is explicitly broken and the inhomogeneity is 
taken to be large, t' <C t. While various sorts of in- 
homogeneities have been found to be widespread in the 
cuprates, it is currently unclear if they are relevant for 
the mechanism of the high temperature superconductiv- 
ity. In part, this is a theoretical issue that hinges on the 
still unsettled issue of what is the superconducting T c (if 



any) of the uniform Hubbard model. Clearly, when there 
is a superconducting groundstate in the highly inhomo- 
geneous limit (t'/t 1), T c is generally an increasing 
function of t' . If T c is small or in the homogeneous 
Hubbard model, there must be an intermediate value of 
t' ("optimal inhomogeneity") at which T c is maximized, 
as suggested in Ref. iJl^Ia- If this is the case, it is sug- 
gestive that the self-organized inhomogeneities found in 
the cuprates may be important for the mechanism of su- 
perconductivity. 
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APPENDIX A: A DERIVATION OF THE 

EFFECTIVE THEORY FOR THE 
CHECKERBOARD HUBBARD MODEL 

Taking the checkerboard Hubbard model as an exam- 
ple, here we would like to show how the effective Hamil- 
tonian can be derived by using (degenerate) perturbation 
theory. To begin with, it is necessary to determine the 
low-energy, reduced Hilbert space. All eigenstates of a 
4-site Hubbard model can be calculated analytically by 
taking advantage of the symmetries of the model. Al- 
though conceptually simple, there are totally 256 eigen- 
states, making the actual calculations involved. Fortu- 
nately, the results have been published by R. Schumann 16 
and explicit eigenstates can be found on his website. For 
U ~ U c , the chosen lowest energy states in each charge 
Q sector are summerized (using Schumann's notation) in 
Table I. 

TABLE I: The eigenstates and eigenvalues of a 4-site Hub- 
bard model in the low-energy Hilbert space. The correspond- 
ing charge secter (Q), total spin (S), spin ^-component (S z ), 
symmetry under tt/2 rotation, and eigenstate number used 
in Schumann's paper are also given, (cosa = [(36t 2 U — 
(7 3 )/27]/[(48i 2 + U 2 )/9] 3/2 , cos/3 = 4i 2 (7/[(16i 2 + U 2 )/3] 3/2 ) 



Q 


E(Q) 


s 


s z 


symmetry 


eigenstate 
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d-wave 


#m 












#46, #70 
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U- \jz2t 2 + U 2 +4 64t 4 +3i 2 U 2 


1 


±§ 


Px ± ip y 
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2 


#50, #74 


2 


U-2y / 48t 2 + U 2 cos(f ) 
3 








s-wave 


#22 



6 



For t' = and 1 ^> x > 0, the groundstate of the 
model is highly degenerate since the doped holes can be 
distributed among the decoupled squares in many ways. 
To an the effective model in powers of t', we employ 
standard degenerate perturbation theory Write Eq. (fTJ) as 
H = Ho + H' (H' represents the t' term) and let V be 
the projection operator onto the subspace spanned by the 
direct product of the states shown in Table I. Then 



H eff =VHV + VH'[1-V]— —[l-V]H'V + .... 

( A1 ) 

For our present study, we keep only terms to first order 
in t'. Specifically, labelling the states in accord with the 
notation of Ref. 16, the various matrix elements of H' 
between the unperturbed groundstates can be expressed 
as follows: the only non- vanishing matrix elements are 
1) elements related to the effective hopping amplitude r, 



*T = <*lllO')| (* 46 (j)l H' |*m0')) l*5oO")) 
= - <*ni(j')| (*bo(j)| H' |* m 0)) |* 46 (i')) 

= (*lll0")l (*70(7')l H ' 1*111 0')) 1*740")) 
= - (*lll0")| (*740)l H' |*m0)) l*7oO")) > 

(A2) 



with j, j' representing nearest-neighbor unit cells; 2) ele- 
ments related to the boson- fermion coupling g, 



9 = (*11iO")I(*m(7')I^' 1*460)) 1*740")) 

= (*2 2 0")l (*lll0')l H' |*460')) 1*740")) 

= (*lll(/)| (*220')l H' |*500)) I*700")) 

= (*220")l (*Hl0')l H' |* 5 o0')) l*700")) 

= - (*lll0")l (*220)l H' |* 74 0)) 1*460')) 

= - (*220")l (*Hl(i)l H' |* 7 40)) 1*460")) 

= - (*lllO")l (*220)l H' |*700)) l*500")) 

= - (*220")l (*mO)l H' |*700)) I*500")) , 

(A3) 



Explicit expressions for these matrix elements as a func- 
tion of U can be obtained using the explicit wavefunc- 
tions given in Ref{H (See Fig. 5). When expressed in 
second quantized form, the resulting effective Hamilto- 
nian is given in Eq. ((3]), above. 
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FIG. 5: The effective parameters, r/t' , g/t' , and A p /t as a func- 
tion of U /t. 
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